Phenomenological theory of the underdoped phase of a high-T c superconductor. 



o 
o 

(N 

X> 
<D 
IX, 

o 

(N 
O 

o 

Oh' 

3 ; 

c3 



^3 

G 

O 

o 



> 
oo 
o 

o 

-I— > 
c3 



i 

o 
o 



X 




A. M. Tsvelik 1 and A. V. Chubukov 2 

1 Department of Condensed Matter Physics and Materials Science, 
Brookhaven National Laboratory, Upton, NY 11973-5000, USA and 

2 Department of University of Wisconsin, Madison, WI 53706, USA 

(Dated: February 3, 2008) 

We model the Fermi surface of the cuprates by one-dimensional nested parts near (0, tt) and (tt, 0) 
and unnested parts near the zone diagonals. Fermions in the nested regions form ID spin liquids, 
and develop spectral gaps below some ~ T* , but superconducting order is prevented by ID phase 
fluctuations. We show that the Josephson coupling between order parameters at (0, tt) and (tt, 0) 
locks their relative phase at a crossover scale T** < T* . Below T" , the system response becomes 
two-dimensional, and the system displays Nernst effect. The remaining total phase gets locked at 
T c < T** , at which the system develops a (quasi-) long-range superconducting order. 
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The copper oxide materials in their underdoped phase 
combine features typical for systems of different dimen- 
sionality. The existence of stripes 0, Q , the absence of 
quasiparticle peaks above T c , the presence of the spectral 
gap near near (ir,0) and (0, tt) at T > T c , and the fill- 
ing of the gap by thermal fluctuations (instead of closing 
of the gap) [3, |4[, are hallmarks of ID strongly corre- 
lated systems. On the other hand, the <i-wave form of 
the pairing gap and existence of curvedportions of the 
Fermi surface (FS) near zone diagonals [5(] are essentially 
two-dimensional features responsible for the physics of 
superconducting cuprates at the lowest energies (see [f|). 
An intriguing aspect of the cuprate physics is the exis- 
tence of the region, above T c , where the system displays 
a strong diamagnetic response 0- This region extends 
to some T** which is smaller than the pseudogap tem- 
perature J 1 *, and has a different doping dependence. 

In this letter we propose a semi-phenomenological 
model which proposes an interpretation of the origin of 
the difference between T c , T** and T* . We assume that 
fermions in antinodal regions near (0,7r) and (7r,0) in- 
teract attractively and form quasi-one-dimensional spin 
liquids (SL) with a large spin gap. These electrons then 
generate low energy collective modes whose dispersion is 
highly anisotropic in space. This assumption is consistent 
with with ARPES experiments [f| showing incoherent 
spectral functions near the antinodes. The assumption 
of strong anisotropy of collective modes is also consistent 
with the 'checkerboard" -like density modulation seen in 
the STM experiments Q , and with the recently observed 
softening of the optical phonons at wave- vectors (4fc F , 0) 
and (0,4fc F ) (Ref. [|). 

The effective interaction in spin liquids is attractive 
and scales to strong coupling. At the other hand, it is 
reasonable to suggest that the coupling at the unnested 
parts and between the unnested and nested parts remains 
weak and repulsive. The increase in doping leads to a 
decrease in the anisotropy which explains opposite trends 
in doping dependence of T* and T**. Without discussing 



FIG. 1: A sketch of the Fermi surface with pairing. Spectral 
gaps open on the nested (blue) regions. 



any microscopic details of the ID SL, we just exploit the 
known fact that interactions there open up quasiparticle 
gaps (pseudogaps) below some crossover temperature T*, 
which is the largest energy scale in our consideration. In 
3D, the gap opening would imply spin or charge ordering, 
in ID the ordering is prevented by strong quantum phase 
fluctuations. 

We consider the Josephson-type coupling between 
the two fluctuating superconducting order parameters 
formed near (0, 7r) and (tt, 0) and demonstrate that be- 
low a crossover temperature T** < T*, fluctuations of 
the relative phase of the two order parameters at (0, tt) 
and (7r,0) acquire a finite mass m ~ j 1 **^ j e ^ fa e re i_ 
ative phase gets locked. This implies that at T** the 
system crosses over into a 2D fluctuational regime with 
enhanced diamagnetic response, and begins displaying 
Nernst effect. We further show that due to topological 
excitations (vortices) in the 2D regime, the (quasi)-long- 
range superconducting order develops only at a smaller 
T c < T** , where the system undergoes XY-model type 
phase transition. Below T c the c?-wave-like spectral gap 
opens on the entire FS due to the proximity effect. 

Our approach bears some similarities to earlier theo- 
ries, but there are important differences. The "stripe" 
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scenario also departs from ID physics, but it assumes 
the existence of two different ID subsystems, separated 
in real space, one with a spin gap, and one without. 
In this situation, the the dominant interaction is within 
each subsystem. In our scenario, the ID SLs near (0, n) 
and (tt, 0) are identical, both having a spin gap, and the 
dominant interaction is between these two SLs. Besides, 
in stripe scenario, the nodal quasiparticles are produced 
by the interaction between the stripes, while we obtain 
nodal quasiparticles due to the interaction between ID 
and curved parts of the FS. In the RVB scenario [loj. 11 1 



SL is supposed to be two-dimensional, and strong phase 
fluctuations at T — emerge from the fact that the 
(Hubbard) interaction is assumed to be larger than the 
fermionic bandwidth unlike in our model where strong 
phase fluctuations are due to quasi- ID behavior. In 
"quantum-critical" scenarios H, 13| , the strongest inter- 
action, responsible for T* , is also assumed to be within 
the antinodal region, and the d— wave superconducting 
gap emerges due to the residual interaction between the 
antinodal regions and the rest of the FS. However, in 
these scenarios, the FS is assumed to be almost circular, 
and the quasi-lD behavior is absent. 

In our model, the FS, shown in Fig. 1, has nested 
regions, and the system can be separated into two quasi- 
lD systems represented by the electronic states close to 
the antinodal regions, and a 2D system represented by 
the states close to the unnested part of the Fermi sur- 
face. The low energy physics is dominated by two gapless 
collective modes originating from critical fluctuations of 
the charge in ID regions, and by nodal quasiparticles. 
Each critical mode is very anisotropic dispersing either 
in the x or y direction. As is always the case in ID, 
such modes describe simultaneously fluctuations of the 
charge density wave (CDW) Acdw and the supercon- 
ducting (SC) Asc — A order parameter . The corre- 
sponding phases Q x and X are related to each other by 
duality d^X = e^d^Ox, (n,v = vt,x), and the ID 
effective action can be equally well expressed in terms 
of X or Q x . Below we use A x ,y as the basic variables 
(subindex refers to fluctuations near (0,7r) and (tt, 0), 
respectively), and also switch from A x ,y to the phases 
X,Y using A x = \A\e lx , A v = \A\e 1Y . At energies 
E « \A\ ~ T*, the order parameter amplitudes are 
frozen, but the phases fluctuate. We assume that the 
most relevant interaction between fermions near (0, tt) 
and (tt, 0) is the Josephson coupling between the phases 
of the two SC order parameters and A y (note that 
the CDW phases & x . y cannot be coupled due to the mo- 
mentum conservation). 

The dynamics of the order parameters X and Y is 
described by the effective action in the form of coupled 



ID XY models: S = j 1/T dr j ^-C u where d is 



~T2 cos [ r O; V) - Y ( x , V + x )} + J cos[X(x, y) - Y(x, y)\ 

Here 2-7r/A is the size of the nested area of the Fermi sur- 
face along a particular direction (experimentally A ~ 4a, 
where a is the lattice constant). For later convenience 
we discretized the gradient term in the action replacing 
it by the cosine. The Coulomb interaction can be ab- 
sorbed into the model parameters as it is screened by the 
ungapped nodal quasiparticles. The last term in §Q de- 
scribes the Josephson coupling between the two SC order 
parameters. We assume that this coupling is repulsive, 
i.e., J > 0, in which case it favors X — Y w tt. 

At J = the action ([1]) describes critical and purely ID 
fluctuations of the phases X, Y. Parameters of the model 
can be expressed in terms of collective mode velocity v 
and the scaling dimension d of the SC order parameter 



(AiTdv) l , t = v/ATrd 



(2) 



which can be extracted from the order parameter corre- 
lation function, e.g., Xxx = (( e ix ( T ^)e- ix (°< ^)) . 

d 

izr./ )- 

\xX 



(3) 



sin (ttTt) + sinh (ttTx/v) 

(Xyy is obtained by interchanging of x and y). 

Deviations from one-dimensionality near (0, tt) or 
(tt, 0) are taken into account by the terms of the form 

C 2 ~ b$(cos[X(x, y) - X(x, y + A)] + 

cos[Y(x, y) - Y(x + A, y)}) + (X, Y -> 4TrQ x . y d) (4) 

These extra terms by itself lead to a finite temperature 
phase transition to a superconducting or a CDW state 
(depending on what type of coupling is more relevant), 
even at J = Q 0. We assume, however, that the cor- 
responding transition temperature is the smallest energy 
scale in the problem, and neglect it, focusing instead on 
the Josephson interaction. 

The coupling of the SC order parameters to the 
unnested part of the Fermi surface is given by 

[j(k x )A x (q) + l(ky)A y (q)]^ k+q/2 yJ+ +q/2 + H.c. (5) 

where 7(...) are phenomenological parameters. This in- 
teraction opens the SC gap on the curved parts of FS 
through the proximity effect when the phases of and 
A y lock. At k x — k y , the couplings 7 are equal, and 
for X — Y = tt (i.e., A x = — A y ) the interaction van- 
ishes. Setting j(k x ) — "fcosk x , we then reproduce the 
d-wave order parameter along the whole FS. When the 
phases X, Y are not locked, the interaction ^ is un- 
able to generate a gap at the unnested part of FS, but it 
can destroy the quasiparticle coherence along the whole 
Fermi surface. Indeed, evaluating the self energy of these 
electrons from (|5l3p in the leading order in 7 we obtained 
that it is local and behaves as 
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It can be shown that the higher order terms are less sin- 
gular in u>. If d < 1 this self-energy exceeds lo 2 coming 
from other, non-singular interactions. Experimentally, 
£(a;) scales nearly linearly with uj above T c . In our phe- 
nomcnological theory, this implies that d « 1/2. 

We now return to effective action (fTJ). Assuming that 
X, Y(x,y) are slowly varying functions of x and y re- 
spectively and replacing t/X 2 terms in ([TJ by quadratic 
functions, we obtain the action quadratic in = 
(X + Y)j\[2. After integrating out we obtain 

the action for <I>~ = (X — Y)/*j2 in the form S e ff — 
J (2#* c eff(k,uj) where 

C eff (k, u>) = [Go ^w, M + Go 1 (w, + 



J(27r) 3 5(k)(5(w) cos(v / 2* 



(7) 



where Go(uj,k) = 2Xd[(ui 2 /v + vk 2 )]^ 1 and means 
Fourier-transform. 

Models of this kind have been discussed before in the 
context of sliding Luttinger liquid phases It was 

assumed there that the theory is renormalizable. By in- 
specting the perturbation theory series in g we have found 
that this is not the case. In this situation, arguments 
based on a supposed relevance of various operators, used 
in 1J] are questionable. In our analysis, we restrict our- 



selves to the summation of the most diverging RPA dia- 
grams. The RPA series for the pairing susceptibility are 
given by 



where i, j = x,y, D = (pJ 2 +tk 2 )(puj 2 +tk 2 : )+m 2 (2pu) 2 + 
tk 2 ) Its poles at iui = E± yield the excitation spectrum: 



£± = y + (™/2) 2 ± ^/ (ttt,/2) 4 + |fc 4 cos 2 (20) (12) 

where k x — k COS 6 hy — k sin 9. The mode E- is gapless, 
as it should be in the ordered state. A substitution of (fTTj) 
into (flO|) yields the self-consistent equation for m solving 
which we reproduce ((9]). 

As we shall see, the vortices transform T** into a 
crossover scale at which fluctuations of <& - develop a gap, 
and the system response becomes 2D. The static diamag- 
netic susceptibility, however, still remains finite, i.e., the 
system does not become a superconductor. To demon- 
strate this, we use the Villain approximation replacing 
all cosine potentials in the effective action {T]) by periodic 
quadratic functions, e.g., 

- cos[X(ar, y) - X(x + X, y)] -> (13) 

- [X(x, y) - X(x + X, y) - 2ttJVi (x, x + A; y)} 2 

We further represent Ni(x,x + X; y) = Ni (x, y) — N\ (x + 
A, y) + Ni(x, x + A; y). The integers Ai j2 live on links of 
the lattice, while N 12 live on its sites. One can easily 
verify that the actions for A\ i2 and Aq j2 decouple from 
each other. Because of the space limitations, we derive 
only the action for N\^- These fields enter S e ff in com- 
bination Ni — N 2 = n. The Lagrangian density is 



Xxx(u; k) = {\xo(u, kc)] 1 - J 2 Xo(u, k y )Y 



(8) 



where Xo(w, k) is the Fourier transform of ([3]). Evaluating 
Xo{k, w) (see Ref. [lj|) and substituting into © we find 
that Xxx(0, 0) diverges at the temperature 



c = P [{xf + on 2 ] + t(d x x) 2 + t(d y Y) 2 + 

m 2 [X(x, y) - Y(x, y) - 2nn(x, y)} 2 



(14) 



Integrating over X, Y fields we obtain the effective action 
for n: 



2tt V A r 2 (l-d/2) 



(9) 



S[n] = 2Tr 2 m 2 T ^ n{-u>, -k)K(oj, k)n(u, k) (15) 



where A = 2itv/X. The divergence implies, that within 
RPA, T** is the true superconducting ordering tempera- 
ture [16( | . Below T**, phase fluctuations acquire a finite 
mass m(T). The spectrum can be obtained by extend- 
ing RPA analysis to the ordered state, i.e., by replac- 
ing cosine Josephson potential by a quadratic function of 

I" = v^" 3 - 7T as Jcos[V2$ (_) ] ~ -tt ($~) 



= J|(cos[V2$ 



(-)ni = 



Jexp 



Plugging this back into ([7]), we get 

({* i (-U,-k)* j (L,,k))) = 

1 ( puj 2 + tk 2 + m 2 m 2 
D(w, k) V m Pu 2 + tkl + m 2 



(10) 



(11) 



K(cj,k) = 



( P UJ 2 +tk 2 x )(pLU 2 +tkl) 



(pcu 2 + tk 2 )(puj 2 + tk 2 ) + m 2 (2puj 2 + tk 2 ) 



At finite T the most important contribution comes from 
zero frequency where the kernel is 



K(0,k)=tk 2 x k' y /(tk 2 x k' y + m 2 k 2 ) 



(16) 



This form of the kernel suggests that the relevant inte- 
ger variable is q(x,y) = d x d y n(x,y) so that n(x,y) = 
J2j ej9(x — Xj) J2k e 'k®(y ~ 2/*)> wnere e j, e 'k are integers. 
Thus from S[n] we arrive to the classical action for the 
Coulomb gas of "charges" q(j,k) = eje' k : 



O'jrf 

S[q] = ^r^2qrln[\r-r'\/b]q r ,, b> 



T 



(17) 
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This form implies that the actual ordering transition be- 
longs to the universality class of XY model (2D one for a 
single layer or anisotropic 3D one for a system of weakly 
coupled layers). The transition temperature then can be 
estimated as 



it irt 
2l 



7TV 

MX 



A 

16d 



(18) 



Below T c the vortices are bound into pairs, both the rel- 
ative and total phase of the two order parameters are 
locked, and the system displays Meissner effect. The cou- 
pling to the unnested part of the FS then gives rise to 
the emergence of the d— wave order parameter all along 
the FS. By construction, the phases X, Y have to be 
locked first, hence T c must be smaller than T**. Para- 
metrically, this requires J to be quite large: JA > A = 
2ttv/X ~ (n/2)v/a. However, numberwise, for d = 1/2, 
we obtain from © and HHD T* * /T c w 70(JA/A), hence 
T** > T c even if JA < A. 

The analysis of the effective action for integer fields 
^1,2 proceeds in a similar fashion. After the integration 
over X, Y we obtain 



The kernel matrix K is such that at zero frequency we 
have the same Coulomb gas action (fP7|) (with the same 
T c !) for the charges q — d y N\ — d x N 2 . In the pres- 
ence of a magnetic field the link variables are modified to 
N\ -> 2eXA x /c + Ni and N 2 -> 2eXA y /c + N 2 , and the 



charges acquire an extra piece q — > q- 



2eX B 



, where B is 



the magnetic field. This implies that these charges are di- 
rectly coupled to a magnetic field, i.e., between T** and 
T c the system develops a strong diamagnetic response 
and displays Nernst effect. 

The existence of the diamagnetic response in a finite 
range above T c is consistent with the measurements of 
the Nernst effect in the cuprates Q. In agreement with 
our theory the temperature up to which Nernst effect has 
been observed is definitely much smaller than T* . Ex- 
perimentally, the onset temperature for the Nernst effect 
follows the same trend as T c and goes down in strongly 
underdoped cuprates implying that the ratio JA/A is in- 
dependent on doping. Our theory also predicts that at 
T > T c , the interaction between nested and curved parts 
of the FS destroys fermionic coherence all along the Fermi 
surface. 

To summarize, in this paper we proposed a phe- 
nomenological model for the low energy physics of under- 
doped cuprates. We assumed that at some large energy 
scale T*, fermions within the nested regions (see Fig.) 
form ID spin liquids. Strong quantum ID fluctuations 
prevent the superconductivity to develop at T*; these 
fluctuations are supressed only at much lower tempera- 
ture T** where the relative phase of the two order pa- 
rameters made of electronic states near (0, 7r) and (tt, 0) 



gets locked, and the system crosses over into a 2D fluc- 
tuational regime with enhanced diamagnetic response. 
Because of the topological vortex fluctuations the total 
phase of the two order parameters becomes locked only at 
even smaller T c < T** , below which the system develops 
a (quasi-) long-range superconducting order. 
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